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We investigated the problem ol the dynamical collapse of a self-gravitating complex charged 
scalar field in Einstein-Maxwell-dilaton theory with a phantom coupling for the adequate fields 
in the system under consideration. We also considered two simplifications of it, i.e., the separate 
collapses of phantom Maxwell and phantom scalar fields under the influence of Einstein gravity. 
One starts with the regular spacetime and leads the evolution through the formation of the horizons 
and the final singularity. We discuss the structures of spacetimes emerging in the process of the 
dynamical collapse and comment on the role of the considered fields in its course. 

PACS numbers: 04.25.dg, 04.40.-b 

I. INTRODUCTION 

Phantom gravitating fields emerge with a kinetic term with the wrong sign. It turns out that this fact implies that 
they are coupled repulsively to gravity. These fields could form ghost condensation at the quantum level, which in 
turn might lead to the modification of gravity at the infrared limit \A\. On the other hand, the present observations 
of the Universe reveal the fact that it may be dominated by an exotic kind of fields with negative pressure causing 
the acceleration. Type la supernovae (Sn la) provide a precise calibrated standard candle with which one can probe 
the expansion of the Universe on large scales. In the late 1990s the independent surveys for distant Sn la announced 
that the high-rcdshiftcd supernovae of this type appeared about 40 percent fainter or equivalently more distant 
than expected in a flat, matter-dominated Universe [2J. These researches were combined with other concerning the 
composition of the Universe such as cosmic microwave background (CMB) observations [3] or analyzes of barionic 
acoustic oscillations (BAO) [4 J. Altogether they do not exclude the possibility that the exotic matter mentioned above 
consists of the fields with the so-called super-negative pressure, i.e., when p < —p (jp and p stand for the pressure and 
density of matter, respectively) [5]. Such fields are dubbed phantom. 

The presence of this kind of fields which contribute negatively to the total energy, i.e., when the violation of the null 
energy condition takes place, gives rise to the interesting possibilities of new solutions of Einstein equations coupled 
to phantom field. Namely, in Ref. [5] the wormhole solutions with phantom scalar field were studied both analytically 
and numerically. Phantom field has also a direct consequence for black hole physics. In Ref. [7 in theories where 
gravity coupled to Maxwell and dilaton fields in such a way that the kinetic energy for either or both of the fields in 
question was allowed to be negative, black hole solutions were obtained. Among all it was found that regular black 
hole solutions with zero or negative mass could be provided. On the other hand, cold black holes with a multiply 
degenerate event horizon and an infinite horizon area were considered in [8|. The new class of static, spherically 
symmetric solutions in Einstein-Maxwell-dilaton theory with phantom coupling for the dilaton and/or Maxwell field 
leading to black holes with single or multiple horizons were elaborated in Ref. [5]. The static multicenter solutions of 
phantom Einstein-Maxwell-dilaton theory and gravitating cr-models obtained via dimensional reduction of phantom 
Einstein-Maxwell, phantom Kaluza-Klein and phantom Einstein-Maxwell-dilaton-axion theories were discussed in 

ma. 

On the other hand, the gravitational collapse and the no-hair conjecture or its mathematical formulation black hole 
uniqueness theorem, stating that stationary axisymmetric solution of Einstein-Maxwell equations relaxes to Kerr- 
Newman spacetime and is characterized by black hole mass, charge and angular momentum |llj . attracted attention 
for many decades of researches. A lot of progress has been also made during last years to answer the basic question 
underlying the theoretical studies of black hole interiors. Quite different description of the inner black hole singularity 
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emerged [T^]- [2]. A full nonlinear description of the inner structure of black holes was presented in Ref.|15j. Then, 
in [16] the spherically symmetric collapse of massless scalar field was numerically examined. It was shown that the 
field dispersed to infinity or collapsed to a black hole. The aforementioned researches were also broaden to the case 
of the nonlinear evolution of the neutral scalar field on a background spacetime of charged black hole [T7]. Ref.fTS] 
was devoted to the numerical solution of spherically symmetric semi-classical scalar field collapse. The addition of 
the effective energy momentum tensor describing evaporation of black hole enabled to consider the influence of this 
process on a critical phenomena present in classical collapse. In |19j the occurrence of the mass inflation phenomenon 
during a dynamical charged gravitational collapse was shown. As far as mass inflation is concerned, it was shown in 
[T4"] that in a realistic RN black hole the Cauchy horizon was unstable due to the fact that it developed a weak null 
singularity caused by the local mass parameter growing indefinitely large. 

The case of the dynamical collapse of spherically symmetric shell of the charged massless scalar field was elaborated 
in Ref.[20 . An external RN spacetime and the inner spacetime bounded by singularity on the Cauchy horizon were 
obtained from numerical calculations. In order to figure out pair creation process in strong electric field in the 
dynamical collapse of electrically charged massless scalar field, the dynamical formation of evaporating spherically 
symmetric charged black hole emitting Hawking radiation and the behaviour of complex scalar field, gauge field and 
the normalized energy momentum tensor during the process in question, numerical approaches were performed |21j - 
[24] . These researches were supplemented by the studies of the behaviour of Brans-Dicke field [25] during gravitational 
collapse of matter. 

Apart from studying the full gravitational collapse on flat background spacetime, the wide range of simulations of 
an accretion onto an existing black hole was also carried out. The most interesting for us exotic matter accreting onto 
the Schwarzschild [26] or Reissner-Nordstrom black hole [27] was modeled as a free scalar field with the opposite sign 
in its energy momentum tensor. 

Because of the fact that the uniqueness theorem for the low-energy string black holes is quite well established 28 , 
the implication of the superstring theory on the dynamical gravitational collapse of charged scalar field was studied 
in [29] . Numerical studies revealed that the dynamical collapse in question led to the formation of Schwarzschild-like 
spacetimes in case of non-zero coupling between the charged scalar field and the fields emerging from the superstring 
theory. 

As was mentioned above, the presence of phantom fields as well as phantom coupling of the ordinary fields, gives rise 
to new possibilities for solutions of generalized Einstein equations. Motivated by the aforementioned analytical and nu- 
merical researches we shall consider the problem of the full nonlinear dynamical collapse of a complex scalar field when 
gravitational interactions take forms of Fiinstem- phantom Maxwell (EM), Einstein-p/ianfom Maiwe^-dilaton (EMD), 
Einstein-Maxwell-p/ianfom dilaton (EMD), and Einstein-p/ianfom Maxwell- phantom dilaton (EAID) theories. We 
will also consider the evolution in Einstein-p/ianfom dilaton (ED) theory, namely the collapse of phantom scalar field 
in the regime of Einstein gravity. 

Our paper is organized as follows. In Sec. II we derive equations of motion describing the gravitational collapse of a 
shell of charged scalar field in Einstein-Maxwell-dilaton theory with phantom coupling. Sec. Ill is connected with the 
numerical scheme applied in our investigations. In Sec. IV we discuss obtained results, while in Sec.V we present the 
conclusions of our investigations. 



II. PHANTOM EINSTEIN-MAXWELL-DILATON THEORY 



The main aim of our researches will be to study the dynamical collapse of a complex scalar field in the string- 
inspired background with phantom coupling of the fields in the theory. Because of the fact that the coupling constant 
of dilaton field to the complex scalar field is unknown, the considered action is provided in the form written in the 
string frame as follows: 




where <p denotes the dilaton field, a stands for the coupling constant between dilaton and charged scalar field if) 
described by the Lagrangian C of the form as 

£ = ~\ (VaV> + *eA a ^ g a0 (V^V* " ieApiP*) ~ &F^F^. (2) 

£i and £2 are dilaton-gravity and [/(l)-gauge field-gravity coupling constants. In the normal case of Einstein-Maxwell- 
dilaton-complex scalar field system they are equal to 1. On the contrary, in phantom case we have £1 = — 1 and/or 
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6 = — 1. After performing the conformal transformation, the action may be rewritten in the so-called Einstein frame. 
It yields 



d X\J — g 



R - 26 (V0) 2 + e 2a ^C^, r , A, e 2 *g a p) 



(3) 



Taking variations with respect to the adequate fields in the theory in question, we obtain the following equations of 
motion: 



V M [ e^F^ 



a + 1 o 

e 20(a+l) 

46 



(Q+1) (v p ip + ieApip\ fvV ~ ieAV) - ^ae 2Q *F 2 = 0, 
ieip* ( W u ip + ieA u iP) - ieip ( Vip* - ieA v ip* 



V 2 %p + ieAP[2Vpip + ieApip) +ieV s A s ip = 0, 



V 2 V>* - ieA p \2Vpip* - ieA tp*j 



ieV 5 A s %p* = 0, 



T„ 



. F, ip, ip*,A). 



The energy momentum tensor T^^, F, tp, tp* , A) is given by 

Tfiui^, F, -0, tp* ,A) = e 2 * (a+1 %M,4>*,A)+T^(F, 

where T^ v {%p, ip*, A) implies 



T^,r,A) 



+ 



ieMA^ S7 v ip* + A v V M V* )- ieip* [An V u ip + A v 



(4) 

(5) 

(6) 

(7) 
(8) 

(9) 
(10) 



V,r V + V,,!- 1 V, ) +i e 2 Afi A v tp 0* + i C(4>,ip*,A)g IMU . 



On the other hand, the explicit form of the Lagrangian £(0, -A) is provided by the expression 

C(ip,tp*,A) = -^[Vpip + ieA^ [V^ip* - ieA^rp*) . 
The energy momentum tensor for £/(l)-gauge and dilaton fields yields 

T^(F, <P) = 6e 2Q ^ Uf^Fu" - \g^F 2 ^j + 6 Uv^A - 9^cP) 2 ^j . 



(11) 



(12) 



In our research we shall start with the regular initial spacetime at approximately past null infinity, compute the 
formation of the horizons and follow the evolution of the dynamically formed black hole to the central singularity. 
One chooses (2 + 2) -spherically symmetric double null coordinate system. It turned out that any coordinate gauge 
transformation of the form u — > f(u) and v — > g(v) will preserve the null character of the retarded and advanced 
coordinates. The line element written in the aforementioned coordinates is provided by [30] : 



ds 2 = —a(u,v) 2 dudv + r 2 (u,v)dfl 2 . 



(13) 



The assumptions about spherical symmetry of the problem imply that the only non-vanishing components of the 
[/(l)-gauge strength field are F uv and F vu . Moreover, another restriction on the gauge potential can be implemented. 
We can consider A u or A v . It turned out that one can restrict his attention to the only one gauge potential due 
to the fact of the gauge freedom A u —> A u + V u 6*, where 8 — J A v dv. Consequently, we are left with the only one 
component of the gauge field being the function of the retarded and advanced coordinates. 
Defining the quantity like in Rcf . 20 : 



A r 2 

Q = 2 A^_r_ 

rt £ 



(14) 
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enables us to separate D-component of the second order partial differential generalized Einstein-Maxwell equations §5§ 
given by 



r 2 e 20(a + l) 

+ ^ ie ( ip* 1>,v -ipip* v )=Q 



into much simpler first order differential equations for A u 



A -9*-o 
Au ' v 2r 2 ~ 



and for Q, which in turn yields 



\ v + 2a <f> tV Q + ^-e 2 * ( W iV - ipip% ) = 0. 



On the other hand, Eq.Q has the explicit form as follows: 

VuVCu + V'.dV' *u + ie A u [ ipip* - ip*ip,, 



' ,uV,v \ ' ,vV,u 



( Q+1 ) p 2*(a+l) _ 
Hi 



6 a 2 Q 2 e 2 ^ 



£i 4r 3 

An inspection of the relations for the complex scalar field ([6|-([7]) reveals that one gets the following: 

ie Q a 2 

r,ui>,v + r v ip u + rip uv +ie r A u ip v + ie r v A u ip H ip = 0, 

r + r ,vi ) * u + r ^*uv -ier A u ip* - ie r „ A u ip* - 16 — - ip* = 0. 

4r 



(15) 

(16) 
(17) 

0. (18) 

(19) 
(20) 



Furthermore, the adequate components of the Einstein tensor and the stress-energy tensor for the underlying theory 
reveal the set of equations 



Sir 



i>,uip*u + ie A u (ip iP* u - ip* ip^ + e 2 A 2 J) iP* 



Cir + \r e 2 ^ Q+1 ) iP, v iP%, 



■ r. 



4r r 
a 2 a r 



e 2a <t> a 2 Q 2 
4^3 : 

Q 2 e 2 ^ a 2 
47 1 



=20(a+l) 



ip tU ip* + ip* u ip tV +ie A u [ip ip* - ip* ip >v 



(21) 
(22) 
(23) 
(24) 



To proceed further, one introduces the new auxiliary variables written in the following forms: 

a a 

s = V>, P = i>,u, Q = ip,v, P = A U , 
k = <p, x = y = 4> }V . 

Moreover, the additional quantities take the forms as 



(25) 



A = T + fg, 

m = fq + gp, 
v = gx + fy. 



(26) 

(27) 
(28) 



It can be easily seen that instead of taking into account two complex fields ip and ip* , one can introduce two real 
fields. Namely, they satisfy the relations ip = ipi + i ip2 and ip* = ip\ — i ipi- Under these assumptions we arrive at 



s = st+is 2 , P = Pi+ip2, <7 = <7i+i<72, 

fi = M2, A*i = fqi + gpi, ^2 = fq2 + gP2- 



(29) 
(30) 
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Just, the system of the second order partial differential equations can be rewritten in the form of the first order one. 
It can be checked that we finally arrive at the set of the first order differential equations written as follows: 



PI : 
P2 : 
P3 : 
PA : 

P6 
P7 : 
P8 : 

El : 
E2 : 

E3& 

£ 4 (1) 
£ 4 (2) 



(He) 



Ml : 
M2 : 



a. u — a c = 0, 
a, v — a d = 0, 
r,u~f = 0, 
r,v -.9 = 0, 

51, u -pi=0, 

52, u~P2 = 0, 

si,v -qi = 0, 

S2,v - 92 = 0, 
k,u 3E = 0, 

k tV -y = 0, 



9,v ~ 


2 di 


f,v + 


A 




r 




A 


9,u + 




r 




A ^ 








A 


d, u - 





= 2fe(a+l) 



Pl 2 + pi + 2 e (si p 2 - s 2 Pi) + e 2 /3 2 (s x 2 + s 2 ) 
(9i 2 + 9 2 2 ) =0, 



2„2 



6e 



2afc Q a 



4 r 3 

2„2 



2afe Q O 



= 0, 



Si a; y + - < 



+ €ix y + - e 



4 r 3 

1 R 2fe(a+1 

4 

. 1 p 2k(a+l) 

4 



Pi 9i + P2 92 + e P (si 92 - s 2 91) 
Pi 9i + P2 92 + e /3 (si q 2 - s 2 9i) 

g a 2 



6e 



2afc 



2 r 4 

l2„2 



, f 2afe<3 O = q 
« e 2 r 4 ' 



4r 
Qa 2 



0, 



n?2,-u +M2 + er/3gi+esi/3g + esi ^ — = 0, 

4r 



rgi.u + /ii-er/3g 2 -es 2 /3g-es 2 



Qa 2 
4r 



0, 



r"92,u + /i 2 + er/3qi+esi/3g + esi — — = 0, 

4r 



a 



2^2 



1 

a + I 

= 0, 



2fe(a+l) 



= 2fc(a+l) 



Pi 9i + P2 92 + e P (s 1 q 2 - s 2 91) 
Pi 9i +P2 92 + e /? (si 92 -s 2 91) 



£1 4 r J 

4i 4 r 3 
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2 a V Q ~ e 2k e r 2 (sx q 2 - s 2 9i) 
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0. 



= 0, 



(31) 
(32) 

(33) 
(34) 
(35) 
(36) 
(37) 
(38) 
(39) 
(40) 

(41) 
(42) 
(43) 
(44) 
(45) 
(46) 
(47) 
(48) 
(49) 
(50) 
(51) 
(52) 

(53) 
(54) 



Let us introduce the quantity of the physical interest. Namely, we define the mass function provided by the relation 

4 



m(u, v) 



1 



f 



fg 



(55) 



It represents the Hawking mass, i.e., the mass included in a sphere of the radius r(u, v). 

The characteristics of all the evolutions under consideration are listed in Table [T] As was explained in detail in 
[TO] the structures of solutions in EMD as well as in phantom-EMD theories may be regarded as equilibrium states 
resulting from the competition and/or cooperation among gravitational, electromagnetic and dilatonic forces. The 
nature of gravitational force is always attractive. The other two forces are either attractive or repulsive depending 
on whether they are phantom or not. It should be kept in mind that in our case we deal with one more constituent 
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TABLE I: Characteristics of evolutions under consideration. The overlining indicates phantom field. Parameters pk and p s 
denote amplitudes of dilaton and electrically charged scalar fields, respectively. 



Type of evolution 


6 


6 


Pk 


Ps 


e 


a 


EM 




-l 




+ o 






ED 


-1 


— 








— 


EMD 


+1 


-l 








0, -1 


EMD 


-1 


+i 








0, -1 


EWD 


-1 


-l 


#0 


7^0 


7^0 


0, -1 



of the system under consideration, namely the complex scalar field. Its influence on spacetime structures is uniform 
in all considered cases, because we do not allow it to be phantom. We shall comment on its role during the collapse 
after a detailed analysis of spacetime structures. 



III. NUMERICAL COMPUTATIONS 



Unfortunately, a detailed analytical investigations of fully nonlinear gravitational dynamics described by the set of 
the conjugate equations (31 )-( 54 1 are impossible. Therefore, numerical methods ought to be used in order to draw 
conclusions about the structure of spacetime emerging in the process of the dynamical collapse. 

In the present paper we used a numerical algorithm containing an adaptive mesh refinement described in [29] . It was 
adjusted for the current issue by taking the parameters £i and £2 into account. The manner of setting the boundary 
and initial conditions for the considered equations remained unchanged in comparison to Ref . [29] . The one-parameter 
families representing the initial profiles of the field functions used for simulations in the present paper are listed in 
Table II] The (/rj)-family refers to the dilaton field with the family constants c\ = 1.3 and C2 = 0.21. The (/s)-family 
describes the electrically charged scalar field. The family constant Vf = 7.5 and the parameter determining the 
amount of the initial electric charge <5 = § (maximal electric charge) |25j . In both cases the free family parameter is 
denoted by p with a subscript corresponding to the particular field. Such choice of the initial conditions is physically 
representative, because the considered collapse is universal, i.e., its course and results do not depend on the types of 
initial profiles. 

There are no analytical solutions neither to the problem in question nor to any of simplified versions of it (apart 
from the extremely simplified, trivial case of empty spacetime). Thus, checking the correctness of numerical code has 
to be based on indirect methods. 

The three tests used for proving the credibility of the code and justifications for carrying them out, are widely 
discussed in [29] . As was reported in the aforementioned paper, the results obtained for the different integration steps 
on the non-adaptive grid display satisfying agreement of an order of 0.01%. Moreover, the code displays the linear 
convergence as is expected for the applied algorithm and errors decrease with increasing grid density. We reached 
the same conclusions for the modified code for all possible combinations of £1,2 shown in Table |TJ Since the figures 
depicting them are analogous to the ones presented in the previous paper, we decided not to present them here. 

The issue of mass (55 1 and charge (14) conservation was also analyzed. We obtained profiles of m and Q versus u 
along ingoing null rays specific to all types of evolutions under consideration. The values of the advanced time were 
chosen in such a way that the null rays did not intersect any horizons for the small values of the retarded time. The 
obtained profiles appeared to be in a qualitative agreement with the ones presented in [29j . For the evolutions leading 
to the formation of a black hole the mass was conserved within 7% and the electric charge within 4.5% in regions 
apart from the vicinities of the apparent horizons. The fluctuations were significant there due to the outgoing fluxes, 
which appear during the refraction of matter on the potential barriers near the gravitational radius of the collapsing 
shell [20 . Mass and charge profiles for non-singular spacetimes also displayed expected behaviour. The values tended 
towards zero for the reason that evolving matter bounced off the non-singular r = and moved towards future null 
infinity leaving empty spacetime 16] . 
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TABLE II: Initial profiles of field functions. 



Family 


Profile 


(Id) 


Pk ■ v • e V <= 2 ) 


(fs) 


Ps ■ sm 2 (tt^) • f cos (tt^) + icos (tt^ + 5) j 



Next, we examined the simplified versions of the problem in question, namely the non-phantom (£i = £2 = 1) 
collapses of neutral and electrically charged scalar fields leading to Schwarzschild and Reissner-Nordstrom spacetimes, 
respectively. The results are consistent with those published, e.g., in Refs. [lT)l [P9l |2"U] . 

The dynamical Schwarzschild spacetime corresponding to the outcome of the collapse of a neutral scalar field |16j 
in fact results from the evolution of dilaton field under the influence of gravitation (£?Z?-collapse). Its structure for 
pk =0.1 is shown in FigJT^i. There is one apparent horizon visible in the spacetime, namely r >v = 0. For v — > 00 
it settles down to an event horizon, which is situated along u = 1.06. It surrounds the central spacelike singularity 
located along the line r = 0, which becomes singular at u = 1.38 (the peak on the diagram). 

The dynamical Reissner-Nordstrom spacetime emerging from the collapse of an electrically charged scalar field (EM- 
evolution) |20| for p s = 0.6 is shown in Figflja. There are two horizons in the spacetime. The outer apparent horizon, 
r i, = 0, after the dynamical part of evolution, when v — » 00, coincides with an event horizon located along u = 0.84. 
The inner horizon, namely the Cauchy horizon is located at future null infinity, i.e., at v — > 00. The line r = 
becomes central spacelike singularity at u = 3.04. 



IV. RESULTS 



As was mentioned in the introduction, the main aim of our research was to investigate the results of the dynamical 
collapse of an electrically charged scalar field in the presence of dilaton field, while one or both types of fields 
are phantom. All the possible configurations were presented in Table [TJ Because of the fact that the structure of 
spacetimes emerging from the dynamical collapse depends on the unknown coupling constant a [29], we take into 
account a = — 1 and a — 0. These values are particularly interesting because they refer to analytical models discussed 
in the literature. Specifically, a = — 1 refers to the low-energy string theory, while when a vanishes the dynamical 
collapse of an electrically charged scalar field in presence of an uncoupled dilaton field takes place. On the contrary, 
the electric coupling constant does not exert any influence on the evolution, so we assume that it is equal to e = 0.5. 
The results of the numerical computations are depicted on Penrose diagrams, i.e., we plot lines of constant r (u, v) in 
the (v, u)-plane. 

The types of the initial field profiles, which were used in our numerical computations are presented in Table |H| The 
(/^-family refers to dilaton field, while the family (fs) is connected with electrically charged scalar field. The free 
family parameter p denotes an amplitude of the particular field. Namely, p s is responsible for the electrically charged 
scalar field, while f>k for the dilaton field. 

At the beginning we will consider the separate evolutions of an electrically charged scalar field with phantom 
coupling of Maxwell field to gravity (£?M-collapse) and phantom dilaton field (_E_D-collapsc) . Then we shall pay 
attention to the evolutions when more than one field is involved (phantom- EM Z?-collapses). In these situations the 
one-parameter cases will be considered at first. It means that the amplitudes of both collapsing fields will be taken 
to be equal and their common value will be denoted by p. Afterwards, in order to complete the physical picture 
of the results of the examined process, for each of phantom-EMD evolutions we shall discuss the behavior of the 
system in the case when the amplitudes of the considered fields differ. To be precise, we describe the changes in 
the spacetime structures appearing when one of the fields has a constant amplitude and the amplitude of the other 
field varies. The constant amplitudes of the electrically charged scalar field and dilaton field were chosen to be equal 
to p s = 0.6 and pk = 0.1, respectively. Such choice was motivated by the fact that these evolutions in the non- 
phantom case lead to the dynamical Reissner-Nordstrom and Schwarzschild spacetimes as was described above. At 
the same time in the phantom cases, as will be explained in sections IV A and IV B[ they lead to Schwarzschild and 
non-singular spacetimes, which are typical for particular phantom evolutions. In view of the bewildering complexity 
of the considered dynamical collapse it will enable us to generalize about the role of the accompanying field (which 



amplitude varies) in the phenomena under consideration. 

A. Einstein-phantom Maxwell collapse 

To begin with, we shall consider the first of the simplified versions of the theory under consideration. Namely, 
we will discuss spacetime structures emerging in the process of the dynamical collapse of the charged scalar field in 
the case of Einstein-phantom Maxwell theory (.EM-evolution). The representative structure is shown in Figj5j The 
parameter characterizing the initial amplitude of the considered field was chosen as p s — 0.6. It turned out that the 
most striking feature of the resulting spacetime is the fact that despite the presence of an electric charge only one 
apparent horizon is observed. The outcome of the considered dynamical collapse is similar to the one in the case of 
Schwarzschild spacetime, which results from the evolution of a self-gravitating electrically neutral scalar field (FigJT^,). 
The single apparent horizon r v — surrounds the spacelike central singularity and settles along an event horizon for 
v — > oo. 

The event horizon is a null hypersurface situated along u — 0.78. In comparison with the case of the dynamical 
collapse of the electrically charged scalar field in non-phantom Einstein-Maxwell case (£7M-evolution) , observed for 
the same parameter p s = 0.6 and depicted in FigJlJ), the emerging event horizon is situated along the smaller value of 
retarded time. As far as the point, where the line r = becomes singular is concerned (peak on the Penrose diagram), 
one can draw the same conclusion. Namely, in phantom dynamical collapse it corresponds to u = 2.98, whereas in 
non-phantom Einstein-Maxwell evolution the singularity arises slightly later as was described above. 

To conclude, one can observe that in the case of Einstein- phantom Maxwell dynamical collapse the structure of the 
issuing spacetime is Schwarzschild-likc and the black hole forms earlier in terms of the retarded time comparing to the 
non-phantom Einstein-Maxwell evolution. The comparison between EM and EM cases allows us to draw an interest- 
ing conclusion concerning the correlation between the attractive/repulsive nature of forces and spacetime structures. 
As was already mentioned, gravitational force is always attractive and supports the collapse. The electrostatic force 
in .EM-evolution is repulsive and the Reissner-Nordstrom spacetime is the outcome of some kind of a compromise 
between the two forces. In contrast, the electrostatic force in iJM-evolution is attractive just as gravitational force. 
In this case a simpler, Schwarzschild-like spacetime emerges. 

B. Einstein-phantom dilaton collapse 

To proceed further, let us consider the second simplification, i.e., dynamical collapse in the case of Einstein-p/ianfom 
dilaton system (_EZ?-evolution), where dilaton field evolves under the influence of gravitational field. This case simply 
refers to the evolution of phantom electrically neutral scalar field. A vast range of the parameters pk characterizing the 
initial amplitude of the field was examined. The results of the numerical calculations for f>k equal to 0.25, 0.5 and 1 are 
shown in Fig[3j An important notion concerning the evolution in question is the fact that for every considered dilaton 
field amplitude the emerging spacetime is non-singular, that is black holes do not form. The only influence of the 
parameter p^ on the resulting spacetime structure is the fact that the bigger value of pk we take the more significantly 
r-coordinate varies in the vicinity of the center of spacetime (pushing out the lines of r = const, in the vicinity of the 
point u = v = 0). The aforementioned behaviour is in complete contrast to the dynamical evolution in Einstein-scalar 
theory (_B£)-evolution). In this case the dynamical collapse of non-phantom scalar field, for appropriately large values 
of the parameter p^ results in a singular, Schwarzschild-type spacetime, as was shown in Fig{T^,. 

Summing it all up, the influence of the phantom nature of the scalar field is even more extreme in the electrically 
neutral case than in the charged one. In Einstcin-p/ianfom dilaton theory black holes do not appear in comparison 
with Einstein-scalar case, where Schwarzschild-type black holes are obtained. It may be also explained in terms of 
the nature of acting forces. Attractive gravitational force cooperates with an attractive force during _EZ?-collapse and 
competes with repulsive dilatonic force in i?Z?-evolution. It turns out that in the latter case gravity is not able to 
overcome the repulsive nature of phantom field. 

C. Einstein-p/iantom Maxwellr Ailaton collapse 

The first theory concerning the collapse of both considered fields, namely the dilaton and electrically charged scalar 
ones, is Einstein- phantom Maa;w;e/Z-dilaton theory (_EMZ3-evolution) . It should be emphasized that there are no 
differences between spacetime structures resulting from the process for the coupling constant a equal to —1 and 
while values of field amplitudes do not differ. In Fig |||the outcome of an evolution in EM D-theoiy for both considered 
values of the dilatonic coupling constant is presented. The value of parameter p was taken to be equal to 0.06 and 
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0.175. For small values of p the spacetime is non-singular. For its bigger values the Schwarzschild-like black hole 
forms. We achieve a single apparent horizon r „ = surrounding spacelike singularity. When w-coordinate tends to 
infinity the apparent horizon coincides with the event horizon of a newly born black hole. 

This situation is also in agreement with the previous analyzes of the nature of acting forces. In this case all of 
them arc attractive, hence for the sufficiently large value of the amplitudes the simplest singular, Schwarzschild-type 
structure forms. 

In order to proceed to considerations of how various fields influence on the dynamical collapse in question, we kept 
an amplitude of electrically charged scalar field p s constant and equal to 0.6 and varied the amplitude of dilaton field. 
The structures of spacetimes for a = —1,0 and pk equal to 0.01 and 0.09 are shown in Fig|5j Then, we performed 
computations the other way round. Namely, we changed an amplitude of electrically charged scalar field p s while 
keeping the other amplitude pt~ constant and equal to 0.1. In Fig[6]we depicted spacetime structures for the amplitude 
of electrically charged scalar field equal to 0.35 and 0.6 for both considered values of dilatonic coupling constant. In 
all the cases we obtained dynamic Schwarzschild spacetimes. 

For small values of a varying amplitude we observed Schwarzschild-like spacetimes like the one shown in FigJlk, 
that is a central spacelike singularity surrounded by one apparent horizon r v = 0, which after a dynamical part 
of evolution settles along an event horizon for v — !• oo. For larger values of a varying amplitude slightly different 
behaviour was revealed, namely two dynamical stages during the evolution were observed. After achieving a constant 
value of it-coordinate an apparent horizon r, v — changes its position to another value of retarded time, which is 
smaller than the previous one. As v — > oo the horizon again settles along u — const, indicating the position of an 
event horizon of a final black hole. Of course a spacelike singularity is situated in the center of it. Such structure also 
corresponds to Schwarzschild-type. 

An increase of pf. results in earlier (in terms of retarded time) appearance of a point, where the line r = becomes 
smg ular (FigJHJ. This increase is also connected with the fact that an event horizon is situated along a smaller constant 
value of it-coordinate. The decrease of a value of retarded time corresponding to a location of an event horizon is 
more significant than the decrease of u referring to the first singular point at r — 0. Within the set of solutions with 
altering amplitude of phantom counterpart of electrically charged scalar field (Figj6]), the changes of p B do not affect 
the location of a point, where the line r = becomes singular. At the same time the location of the event horizon 
moves towards lower values of it-coor dina te when p s increases. 



The energy momentum tensor (|9j-(12) for EM D-theory in the case of a = —1 has the following form: 

T^(F^) = -e-^T^W+T^ty), (56) 



while for a — 0: 



Tfu/iF, (ft) = -T^(F)+T^{4>), (57) 
T^F,i>^\A) = e^f^(4,^\A)+T^{F,4>), 

where T M „(t/>, ip* , A), T^(F) and T^ v ((ft) denote the parts connected with electrically charged, Maxwell and dilaton 
fields, respectively. The coexisting effects of the dilaton field and phantom counterpart of an electrically charged scalar 
field result in an appearance of Schwarzschild black hole. Even though the changes of the amplitudes of particular 
fields change their relative impact on the collapse, its outcome is in all cases similar. 

It is worth mentioning that analytic studies of singular static configurations in EMD-theory were conducted in 
Refs. pj |S] . It turned out that in the considered case the solutions are Schwarzschild-type, i.e., the spacetimes consist 
of one horizon surrounding central singularity. This conclusion was confirmed by our computations also in the case 
of dynamical collapse. 



D. Einstein-Maxwell-p/taniom dilaton collapse 

The next type of the theory we shall be interested in is Einstein-Maxwell-p/taniom dilaton gravity (FAfD-cvolution) . 
In this theory the electrically charged scalar field coupled to phantom dilaton field collapses gravitationally. The plots 
bounded with the dynamical collapse in FAfU-evolution, for the dilatonic coupling constant equal to —1 and are 
shown in Figs(7]and[8j respectively. In our considerations, for the case of a = — 1, we have chosen the parameter p to 
be equal to 0.5, 0.85 and 1.1. When the coupling constant a is equal to zero p equals to 0.1, 0.8 and 1.1. 

For all values of the family parameter p the non-singular spacetimes emerge for both considered values of a. For 
small values of the parameter p there are no horizons in spacetime. Refining our numerical studies to the larger values 
of p, we obtain an apparent horizon r „ = for small values of retarded time. For a = —1, it settles along constant 
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it-coordinate for v — ¥ oo, while for a — it forms an open loop with ending points at v — > oo. It happens that for 
even bigger values of the amplitudes another type of an apparent horizon, namely r „ = 0, emerges. We observe both 
r u = and r v — types of the apparent horizons. However, their presence in spacetimc is not connected with the 
appearance of the singularity. 

For a = — 1 the results are depicted in FiglT) One observes that the appearance of r :U = and r „ = apparent 
horizons for small values of M-coordinate develops a gap between r — const, lines at v — > oo. The adjacent lines move 
away very quickly. As was mentioned above, for large values of parameter p both types of apparent horizons emerge. 
For the coupling constant a = the results are shown in FigjS] One gets that for large values of the parameter p and 
small values of the u-coordinate as well as for large values of both retarded and advanced times, both r u = and 
r jV = apparent horizons exist in the non-singular spacetime. 

In the considered EM D-theory the repulsive electromagnetic and dilatonic forces compete with the attractive 
gravitational force after all leading to non-singular spacetimes with single or multiple horizons. 

In order to get a deeper insight into the physical processes taking place during E'MD-collapse we shall vary 
amplitudes of phantom dilaton and of electrically charged scalar field. We commence with altering an amplitude of 
phantom dilaton while an amplitude of the other field remains constant and equal to p s — 0.6. For a = — 1 the 
structures formed when pk equals to 0.15, 0.2 and 0.25 are shown in Fig|9] For small values of the parameter pk we 
get the singular spacetime with two kinds of apparent horizons, but when the amplitude of phantom dilaton grows 
the spacetime becomes non-singular. The growing amplitude of phantom dilaton field suppresses emergence of a 
singularity. This result seems to be intuitive because phantom field contributes negatively to the total energy of the 
system. 

As was mentioned above, for small values of the phantom dilaton amplitude the spacetime became singular but with 
two types apparent horizons r u = and r v = in the domain of integration. For values of amplitudes not exceeding 
0.15 there are two branches of r, = horizon. The first one settles along an event horizon for v — > oo. With an 
increasing value of the phantom dilaton field amplitude, the singularity as well as the horizon appears at increasing 
u-coordinate. The second branch of the horizon r :V = lies inside the first one and it begins at the point, where the 
line of singular r = ends. The apparent horizon r u — is situated inside the second branch of r v = horizon 
and begins at the same point. For small values of parameter pk these two lines almost coincide. The structure of an 
area beyond the horizon r M = is impossible to interpret fully and clearly, because calculations cannot be carried 
out until r = and they end at some positive value of r. Such structure with two branches of r iV = apparent 
horizon and an apparent horizon r u — inside them may be interpreted as a dynamical wormhole [24 . On the 
other hand, for greater values of the parameter f>k (larger than 0.2) the two branches of the horizon r, = form 
a loop with two ends coinciding with ending points of singular part of line r = 0. The horizon r „ = remains 
outside it and continuously tends towards smaller values of u-coordinate as v — > oo. This structure corresponds to an 
unphysical situation of naked singularity in spacetime [24j . For even bigger values of amplitudes of dilaton field, from 
approximately pk — 0.25, the spacetime is non-singular wit h no visible horizons. 



A close inspection of the energy momentum tensor (|9j)-( 12 1 for EM D-theory in the considered case of a = —1 reveals 
that one gets the following: 

= e-^T^CF)-!^ (</>), (58) 

It can be seen that the factor e~ 2 ^ diminishes the influence of the Maxwell field on the collapse, i.e., the bigger 
amplitude of phantom dilaton field we consider the smaller influence is of the £/(l)-gauge field we obtain. When 
pk increases the strong repulsive effect of the phantom dilaton field begins to play the dominant role leading to 
non-singular spacetimes. 

Different situation takes place for Einstein-Maxwell-p/wmfom dilaton collapse without dilaton coupling, that is for 
a — 0. The structures of the emerging spacetimes for pk equal to 0.05 and 0.8 are shown in Fig |10| The spacetimes are 
non-singular for all values of the parameter pk ■ It means that the presence of phantom dilaton field, not interacting 
with the rest of matter in spacetime, prevents black hole formation. The spacetime is non-singular with both types 
of horizons at first and then, for pk exceeding 0.7, only one open loop-shaped apparent horizon r u = is observable 
for the small values of the w-coordinate. 

On the other hand, the energy momentum tensor for a — is provided by 

V(F,0) = T M1/ (F)-T^ (<£), (59) 



Eq. ( 59 1 reveals that the growth of the phantom dilaton field amplitude plays the dominant role while the factor e 2 ^ 
also enlarges the influence of the electrically charged scalar field. Therefore due to these facts the domination of the 
phantom dilaton field seems to be a key factor in the process of the dynamical collapse for this case. 
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The next set of solutions referring to Einstein-Maxwell-p/wmfom dilaton theory involves spacetimes emerging from 
a collapse when an amplitude of phantom dilaton field is constant and an amplitude of electrically charged scalar field 
varies. The value of the constant amplitude was set as pj. = 0.1. In Fig{TT]the structures of spacetimes emerging 
from the collapse for a = — 1 were depicted. The amplitudes of the scalar field p s were chosen to be equal to 0.3 and 

0. 6, respectively. One can observe that for small values of the electrically charged scalar field amplitude the repulsive 
character of phantom dilaton field combined with small amount of the electrically charged scalar field results in the 
nonsingularity of the emerging spacetime. On the other hand, when the role of phantom dilaton field is diminished, 

1. e., the amplitude of electrically charged scalar field grows, one attains the singular spacetime with two branches of 
r, v = apparent horizon - the outer and the inner one as well as an apparent horizon r u = located within them. 
One remarks that a dynamical wormhole structure appears. Namely, the emerging black hole gains phantom dilaton 
matter violating the null energy condition and this phenomenon causes formation of a wormhole. 

As may be inferred from Eq.(58) the growth of an electrically charged scalar field amplitude is able to neutralize the 
diminishing effect of the factor e~ 2< ^ coexisting with the repulsive tendency of phantom dilaton field. It finally leads 
to the formation of a dynamic wormhole. 

The structures of spacetimes emerging from .EAfD-evolution with constant value of an amplitude pf. and dilatonic 



coupling constant a — are shown in Fig 12 The amplitude of charged scalar field p s was set as equal to 0.2 and 0.7. 
The spacetimes for all values of a varying amplitude are non-singular. For small values of it there are no horizons 
visible in the domain of integration. In turn, for bigger values of p s apparent horizons r v = and r u = appear in 
spacetime, but none of them reach the line r = 0, which for this reason does not become singular. The emergence of 
the apparent horizons is accompanied by squeezing the lines of r — const, in the vicinity of the line r = 0, in the area 
lying abo ve t he r „ = horizon in terms of the u-coordinate. 

From Eq.(59) we can observe that the factor containing phantom dilaton field e 2 ^ plays an important role together 
with the growth of the electrically charged scalar field amplitude in the domain of integration. This is closely related 
to the fact that phantom dilaton field suppresses appearance of black hole, while an increasing amplitude of the 
electrically charged scalar field relatively to dilaton field amplitude forces the formation of horizons in the underlying 
spacetime. 



E. Einstein-phantom Maxwellrphantom dilaton collapse 



The last examined case is Einstein-p/iantom Maxwell-phantom dilaton theory (EM _D-evolution) . The resulting 
spacetimes formed during the process of the dynamical collapse are depicted in Figs |13| and |14| The plots were made 
for the dilaton coupling constant equal to —1 and 0, respectively. For a = — 1 the amplitudes of the considered fields 
were chosen to be equal to 0.3, 0.4 and 0.5. In the case of a = 0, the parameter p was equal to 0.2 and 0.5. 

For both values of a the spacetimes for small values of p are non-singular without any horizons. Similarly to the 
previously analyzed dynamical collapse, in the majority of cases both types of apparent horizons, namely r v = and 
r n = are visible in the domain of integration. Moreover, the appearance of the apparent horizon r „ = for larger 
values of the free family parameter always ensues from the formation of singularity in the considered spacetime. What 
is more, the bigger value of the parameter one considers the earlier in terms of retarded time the singularity at r = 
appears. 



In Fig 13 plotted for the coupling constant a = — 1, one observes a loop-shaped part of the apparent horizon r iV = 
near the peak of the line r — 0, where it becomes singular. The size of a loop-shaped apparent horizon decreases with 
increasing field parameter p. What is more, the loop surrounds the singular part of the r = line. Outside the loop 
the singularity is absent or impossible to localize. When ^-coordinate tends to infinity, the apparent horizon r„ = 
settles along constant it-coordinate. Moreover, the other apparent horizon begins at the point, where the singular 
r = ends and tends towards smaller values of retarded time as v —> oo. As was already stated in the previous 
section, such structure corresponds to naked singularity in spacetime. 



In Fig 14 we presented the results for the uncoupled case, that is for a = 0. One notices that for the large values of 
parameter p the spacetime relaxes to the one containing naked singularity. There is a loop-shaped apparent horizon 
surrounding the central spacelike singularity and it tends towards u = const, as v — > oo. Additionally, there are two 
branches of the horizon r iU — 0. One of them is connected with the naked singularity. The other one forms an open 
loop similar to the ones observed in spacetimes emerging from the iSMD-evolutions conducted for a = 0. 

In conclusion it may be stated that spacetimes emerging from the collapse in Einstein-p/ianfom Maxwell-phantom 
dilaton theory with equal amplitudes of both fields are either non-singular or they contain naked singularities. In this 
case only the dilatonic force is repulsive, the other two forces have attractive character. 

Let us turn our attention to the evolutions with a constant amplitude of phantom electrically charged scalar field 
equal to p s = 0.6. For a — —1 we considered the values of phantom dilaton field amplitude pk equal to 0.1 and 
0.3. The spacetimes emerging from these evolutions are presented in Fig[l5j For all values of pk the structure of the 
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spacetime resembles Schwarzschild one, i.e., there is an apparent horizon r .„ = coinciding with an event horizon at 
v — > oo and surrounding central spacelike singularity. The peak, where the line r — becomes singular practically 
does not change its it-location as pk increases. On the contrary, the value of it-coordinate corresponding to the location 
of the event horizon increases while the amplitude of phantom dilat on fi eld becomes bigger. 



In the considered case of a = — 1 the energy momentum tensor (J9|)-(12) for EM D-theory is given by: 

T„ v {F,<t>) = -e-^T^W-T^ty), (60) 



Contrary to the case described by relations (58 1, where we observe dynamic wormhole spacetimes, now the phantom 
coupling of Maxwell field causes that the outcome is similar to the Schwarzschild spacetime. It should be emphasized 
that it happens despite of the influence of phantom dilaton field, which previously prevented black hole formation due 
to its repulsive nature. 

On the other hand, Fig |16| illustrates spacetime structures resulting from a collapse in Yimstcm- phantom Maxwell- 
phantom dilaton theory without dilaton coupling, i.e., for a — 0. The values of phantom dilaton field were set as 
equal to 0.05 and 0.2. For small values of ph the spacetime is singular with one apparent horizon r „ = surrounding 
singularity visible in the domain of integration. Nearby the point, where the line r = becomes singular the apparent 
horizon almost coincides with the singularity. At some point it moves away from the line r = towards slightly 
smaller values of ^-coordinate and then settles along constant it for v — > oo. Since it becomes null hypersurface there 
we may predict that its location at large values of advanced time indicates the location of the event horizon. Although 
the structure in the vicinity of the point, where the line r = becomes singular (beyond the event horizon) is different 
from the one presented in Fig[T^, we conclude that it may be interpreted as Schwarzschild-like. The emerging object 
consists of a central spacelike singularity surrounded by one apparent horizon coinciding with an event horizon at 
v — > oo. The length of a part, where the apparent horizon is very close to r = decreases with an increase of an 
amplitude of phantom dilaton field. The location of the point, where the line r = becomes singular slightly decreases 
in terms of u while the varying amplitude increases. For larger values of a varying amplitude, naked singularity and 
an additional open loop-shaped branch of an apparent horizon r tU = appear in the spacetime. 
Consider now the energy momentum tensor for the studied a — case. It implies 

T^{F,4>) = -T^{F)-T^{4>), (61) 
T^(0,F,i;,r,A) = e 2 *f^,r,A)+T^(F,cj)). 



Having in mind Eq.(59), the case when no formation of a black hole was indicated, now we observe that phantom 
coupling of Maxwell field provides the singularity in the emerging spacetime and the formation of a Schwarzschild-like 
black hole or naked singulatrity. Similarly to the case of a = —1 black hole formation takes place in spite of the 
presence of the repulsive phantom dilaton field. 



The next set of solutions corresponding to i?M_D-theory refers to the computations conducted for a constant 
phantom dilaton field amplitude p^ =0.1 and altering amplitude of phantom counterpart of an electrically charged 
scalar field. The spacetimes ensuing from the collapse for dilatonic coupling constant a = — 1 and values of an 



electrically charged scalar field p s equal to 0.25 and 0.5 are shown in Fig 17 For small values of an amplitude p. 



the spacetime is non-singular, while for larger values it becomes Schwarzschild-type. There is one apparent horizon 
r „ = indicating the location of an event horizon for v — > oo. It surrounds central spacelike singularity like in FiglTk. 
The damping effect of e -2 ^ in (60) is dominant at first and the phantom dilaton field plays the most important role. 



But with the passage of the increase of the phantom counterpart of the electrically charged scalar field amplitude 
the situation reverses and this field manages to force black hole formation. Just one can conclude that instead of 
the diminishing factor e^ 2 ^, the sign of energy momentum tensor of [/(l)-gauge field combined with the increasing 
value of the amplitude of the electrically charged scalar field are of the crucial importance for the formation of 
Schwarzschild-like black hole. 

The diagrams of spacetimes for a = are illustrated in Fig jl8| The amplitude of the phantom electrically charged 
scalar field was equaled to 0.35 and 0.65. For small values of p s the spacetime is non-singular. The structures obtained 
when the amplitude p s increases are Schwarzschild-like with an additional branch of an r tU = apparent horizon at 
small values of the retarded time. 

On the other hand, a close inspection of (61) reveals that together with the dominating factor e 2 ^, the opposite sign 



of the Maxwell term again occurs decisive for the singularity to appear in the domain under inspection. 

As was noticed in the above subsection, the very familiar situation as described in EM D- theory, i.e., black hole 
solution with effectively one event horizon in the case of a 2 — > 1 [7 , takes place in EMD-gr&vity. This analytical 
result obtained in the static case was confirmed in our simulations of the dynamical collapse. 
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V. CONCLUSIONS 



In our paper we studied the problem of the dynamical collapse of charged scalar field in the Einstein-Maxwell-dilaton 
theory with phantom coupling of the adequate fields. These fields are allowed to violate the null energy condition. 
The landscape of static and spherically symmetric solutions for the system under consideration abounds in eleven 
different not only asymptotically flat spacetimes related to sixteen various causal structures [5j . Such a rich collection 
of solutions aroused our interest especially in the context of the dynamical evolution. 

At first, it was revealed that in Einstein- phantom Maxwell evolution the structure of the resulting spacetime re- 
sembles Schwarzschild-like spacetime despite of the presence of the electrically charged [/(l)-gauge field. In this 
case phantom nature of the electrically charged scalar field forces Schwarzschild-like structure formation instead of 
Reissner-Nordstrom-like in the non-phantom case. What is more, the resulting event horizon in the case of EM 
dynamical collapse emerges slightly earlier, in terms of the retarded time, than in the case of the dynamical collapse 
in .EM-theory. Such structure remains in agreement with the one described in non-dynamical case pj [9] ■ 

Then, taking into account Einstein- phantom dilaton theory, one observes that for every considered amplitude of 
the dilaton field the resulting spacetime is non-singular. There is no formation of a black hole. Such behaviour 
is in complete contrast to the dynamical collapse in Einstcin-dilaton (Einstein-scalar) evolution, which leads to the 
formation of Schwarzchild-like spacetime 16 j. In this case the phantom nature of scalar field prevents a singularity 
as well as event horizons formation due to the repulsive character of the field. 

On its turn, in Einstein- phantom Maxwell-dilaion case there are no qualitative differences among the structures of 
the emerging spacetimes for the different values of the coupling constant a. One arrives at the Schwarzschild-like 
black hole formation, with a single horizon surrounding spacelike singularity. Dilaton field and phantom counterpart 
of electrically charged scalar field cooperate in Schwarzschild-type structure formation. These dynamically obtained 
structures only partially confirm the static solutions for EM D-theory obtained in [9]. From a wide variety of pos- 
sibilities during the dynamical evolution the simplest one with only one horizon surrounding spacelike singularity is 
preferred. 

As far as Einstein-Maxwell-p/iantom dilaton theory is concerned, the obtained structures are far more complicated 
than in the previous cases. First of all, two types of apparent horizons are observed in the spacetime, namely lines of 
r_ u = and r jV = 0. Moreover, the appearance of an apparent horizon is not necessarily connected with the emergence 
of the singularity. For the dilatonic coupling constant equal to —1 dynamical wormhole spacetimes are observed in 
case when an amplitude of phantom dilaton field is considerably smaller than an amplitude of an electrically charged 
scalar field. In fact one has a formation of the black hole which gains phantom dilaton field which violates null 
energy condition. This phenomenon provides a formation of the wormhole structure. For a — the spacetimes are 
non-singular. It means that an uncoupled phantom dilaton field prevents black hole formation due to its repulsive 
character. It seems that in the case of EM D-co\la,pse the diversity of solutions manifests itself most clearly. Obviously, 
as was noticed in [TU] , multiple horizons separating distinct areas are observed. Moreover, as was anticipated in 
static cases [JJ [5] , wormhole structures in the case of phantom scalar field and non-phantom vector field were obtained 
also during the dynamical collapse. Interestingly, our computations also provided an evidence of the existence of the 
dynamical interconversion between black holes and wormholes, which results from the dynamics of the phantom field 
in the considered spacetime [3T] . 

Finally we note that, in the case of Einstein- phantom Maxwell-phantom dilaton theory, also both types of the 
apparent horizons are visible in the domain of integration. The obtained structures correspond to the non-singular 
or Schwarzschild-type spacetimes or they contain naked singularities. Naked singularities mainly appear when the 
amplitudes of both collapsing fields are equal. On the other hand, the dynamical collapse results in Schwarzschild-like 
spacetime when an amplitude of the dilaton field is considerably smaller than the other one. Hence it may be stated 
that the considered i?MD-evolution again leads to the simplest structure with a singularity surrounded by an event 
horizon among those presented in [9]. 

During our numerical simulations of a dynamical collapse in phantom-EMD system we obtained an extensive 
collection of the spacetime structures. They are summarized in Table III together with the structures emerging from 
the collapse without phantom coupling. An analysis of energy-momentum tensors and the nature of forces acting in 
the considered system combined with the study of the emerging spacetime structures, led us to the general conclusions 
concerning the role of particular fields and the dilatonic coupling during the collapse in question. 

It turned out that the non-phantom dilaton field supports the formation of Schwarzschild-type singular spacetimes. 
On the other hand, its phantom counterpart suppresses black hole formation due to its strongly repulsive nature. It also 
facilitates the formation of multiple horizons not necessarily connected with the singularity. The non-phantom Maxwell 
field seems to be less repulsive than phantom dilaton field and thus it favors Reissncr-Nordstrom-type spacetime as 
an outcome of the dynamical collapse. What is more, phantom Maxwell field, similarly to the dilaton field, supports 
the creation of the simplest singular, Schwarzschild-like structure despite its phantom coupling to gravity. 
Considering the role of a complex scalar field during the collapse, it seems to act just as a carrier of the Maxwell field. 
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TABLE III: Spacetime structures emerging from the phantom gravitational collapse. The overlining indicates phantom field. The 
structures of spacetimes resulting from the evolutions in non-phantom gravity are summarized at the bottom for comparison. 
Symbols: A - attractive, TZ - repulsive, ns - non-singular spacetime, S - Schwarzschild-type spacetime, RN - Reissner- 
Nordstrom-type spacetime, WH - wormhole, NS - naked singularity. Superscripts: * - additional horizons in spacetime, 2 ' - 
two-stage collapse. Arrows: f - an increase of a particular amplitude, — > - the direction of changes while the amplitude varies. 
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It does not influence the obtained structures distinctly. Such situation is in agreement with the previously examined 
simpler case of the collapse of electrically charged scalar field in Einstein gravity [20] . 

Considering the dilatonic coupling it may be revealed that a = — 1 diminishes the influence of Maxwell field on 
the evolution and thus enhances the role of the dilaton field in its course. The dilatonic coupling constant equal to 
zero acts the other way round, i.e., it reduces the role of dilaton field during the collapse and indirectly enlarges the 
meaning of Maxwell field by acting on the complex electrically charged scalar field. 
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FIG. 1: (color online). Penrose diagrams of the dynamical (a) Schwarzschild and (b) Reissner-Nordstrom spacetimes obtained 
for the field amplitudes equal to = 0.1 and p a = 0.6, respectively. 




FIG. 2: (color online). Lines of constant r in the (ww)-plane for i?M-evolution with the electric coupling constant e = 0.5. The 
family parameter is set to be p 3 = 0.6. 
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FIG. 4: (color online). Lines of constant r in the the (uu)-plane for _BA/D-evolution with the electric coupling constant e = 0.5 
and the dilatonic coupling constant a equal to —1 and 0. The family parameter is taken to be (a) p = 0.06 and (b) p = 0.175. 
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FIG. 5: (color online). Lines of constant r in (uit)-plane for the EMD evolution with a varying dilaton field amplitude. The 
coupling constants have the same values as in Fig|l] The family parameter for the electrically charged scalar field p a — 0.6, 
while for the dilaton field (a) pk = 0.01 and (b) pk = 0.09. 
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FIG. 6: (color online). Lines of constant r in the («tt)-plane for the EMD evolution with varying electrically charged scalar 
field amplitude. The coupling constants have the same values as in Fig(4] The family parameter for the dilaton field is equal 
to pk — 0.1, while for the electrically charged scalar field field is set to be (a) p s = 0.35 and (b) p s = 0.6, respectively. 




FIG. 7: (color online). Lines of constant r in the (uit)-plane for _EMD-evolution with the electric and the dilatonic coupling 
constants equal to e = 0.5 and a = — 1, respectively. The family parameter is equal to (a) p = 0.5, (b) p — 0.85 and (c) p — 1.1. 
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FIG. 9: (color online). Lines of constant r in («u)-plane for EMD evolution with a varying dilaton field amplitude. The electric 
and dilatonic coupling constants are equal to e = 0.5 and a = —1, respectively. The family parameter for the electrically charged 
scalar field is equal to p s — 0.6, while for the dilaton field is set as (a) pt = 0.15, (b) pk = 0.2 and (c) ph = 0.25. 
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FIG. 11: (color online). Lines of constant r in the (vu)-plane for EMD evolution with a varying electrically charged scalar field 
amplitude. The electric and dilatonic coupling constants are equal to e = 0.5 and a = — 1, respectively. The family parameter 
for the dilaton field is pk = 0.1, while for the electrically charged scalar field field is chosen to be (a) p s = 0.3 and (b) p s = 0.6. 
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FIG. 12: (color online). Lines of constant r in the (nu)-plane for the evolution described in Fig |ll| with the dilatonic coupling 
constant equal to a = and the family parameter for electrically charged scalar field field (a) p~ s — 0.2 and (b) p s = 0.7. 




FIG. 13: (color online). Lines of constant r in the (uit)-plane for I?M D-evolution with the electric and dilatonic coupling 
constants equal to e = 0.5 and a = —1, respectively. The family parameter is set (a) p — 0.3, (b) p — 0.4 and (c) p = 0.5. 
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FIG. 14: (color online). Lines of constant r in the (vit)-plane for iJAfD-evolution with the electric and dilatonic coupling 
constants equal to e = 0.5 and a = 0, respectively. The family p arameter is given by (a) p = 0.2 and (b) p = 0.5. 
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FIG. 15: (color online). Lines of constant r in the (wu)-plane for EMD evolution with a varying dilaton field amplitude. 
The electric and dilatonic coupling constants are equal to e = 0.5 and a = — 1, respectively. The family parameter for the 
electrically charged scalar field p s = 0.6, while for the dilaton field (a) f>k = 0.1 and (b) pt = 0.3. 
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FIG. 16: (color online). Lines of constant r in the (tm)-plane for the evolution described in Fig 15 with the dilatonic coupling 
constant equal to a = and the family parameter for dilaton field (a) pk = 0.05 and (b) pk = 0.2. 
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FIG. 17: (color online). Lines of constant r in the (iiu)-plane for EMD evolution with a varying electrically charged scalar field 
amplitude. The electric and dilatonic coupling constants are equal to e = 0.5 and a = — 1, respectively. The family parameter 
for the dilaton field is chosen as pt = 0.1, while for the electrically charged scalar field field equals to (a) p s — 0.25 and (b) 
p s — 0.5, respectively. 
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FIG. 18: (color online). Lines of constant r in the (tm)-plane for the evolution described in Fig 17 with the dilatonic coupling 
constant equal to a = and the family parameter for the electrically charged scalar field field equal to (a) p s = 0.35 and (b) 
p s = 0.65. 



